The role of microscopic fluctuations in transition predictior£| 
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The commonly accepted description of transition to turbulence in shear flows requires the presence 
of an external source of disturbances that get amplified by an essentially linear mechanism up to the 
point where breakdown to turbulence occurs. Microscopic fluctuations are shown here to provide 
just the right amount of initial disturbances to match the predictions of linear stability theory. 
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I. INTRODUCTION 

Things as different as Brownian motion and the in- 
put noise of an electronic amplifier (and Kirchhoff's law 
of thermal radiation, although not generally mentioned 
in this connection) have their common explanation in the 
fluctuation-dissipation theorem, which stipulates that ev- 
ery macroscopic dissipation phenomenon is accompanied 
by a random forcing induced by the underlying micro- 
scopic mechanics, with a spectrum that is a universal 
function of temperature. In fluid dynamics the appropri- 
ate forcing, in the form of a random stress tensor, was 
added to the Navier-Stokes equations by Landau & Lifs- 
chitz [l| and further discussed by Fox & Uhlenbeck 0]. 

In almost the same years, the theory of linear instabili- 
ties of the Navier-Stokes equations gradually evolved into 
an engineering tool for transition prediction. Being able 
to determine the position of transition from laminar to 
turbulent flow has a large practical impact on the design 
of aircraft and many other man-made artifacts; at the 
same time it is a problem that still defeats a complete 
solution, in spite of universal agreement that the Navier- 
Stokes equations constitute an adequate physical model. 
It is significant that, in an often quoted sentence, Richard 
Feynman :3j chose a transitional flow to produce a didac- 
tic example of the inscrutability of differential equations: 
The next great era of awakening of human intellect may 
well produce a method of understanding the qualitative 
content of equations. Today we cannot. Today we can- 
not see that the water flow equations contain such things 
as the barber pole structure of turbulence that one sees 
between rotating cylinders. 

The starting point of numerical transition prediction 
was the (rather inexpressively named, but so known ever 
since) e^ method of Smith & Gamberoni jj] and Van In- 
gen . The crux of it was the realization that shear flows 
behave as ampliflers (in a more recent language, they 
support convective instabilities), accepting small distur- 
bances of various origins at some upstream location and 
amplifying them downstream to the point where tran- 
sition occurs. The original version of the method only 
computes the linear-ampliflcation factor e^ of ToUmien- 
Schlichting waves, from one of a few alternative quasi- 



parallel approximations of the Navier-Stokes equations, 
and uses the value of the exponent as a transition crite- 
rion (llj . experiments having established that transition 
occurs at the location where w 9-^ 10 under the small- 
est environmental disturbances, or earlier otherwise. 

In an effort to improve upon the e^ method many re- 
searchers have focussed on the receptivity problem, i.e. 
the coupling of external disturbances to the boundary 
layer that precedes the ampliflcation proper. For a re- 
view see, e.g., Saric 0|. Different types of disturbances 
have been studied in a relatively large literature: acoustic 
waves, freestream vorticity, wall roughness and the pres- 
ence of a leading edge are all considered to be sources 
of the boundary-layer instability mode that, once am- 
plified, triggers transition. (Nonlinear mixing of at least 
two of these sources is generally necessary to produce the 
right phase speed.) However, a stumbhng block to the 
application of receptivity theory to transition prediction 
has been the difficulty of specifying realistic amplitudes 
and spectra for all these external disturbances. For this 
reason, most if not all engineering computations are still 
today performed by the e^ method without any recep- 
tivity accounted for. 

Despite microscopic fluctuations having been given at- 
tention by authors who pursue a general theory of tur- 
bulence {e.g. 0) and, recently, in DNS of Couette flow 
Q , surprisingly they have never been considered as can- 
didate disturbance sources in the context of transition 
prediction. Yet the e''^ method views the boundary layer 
as an amplifler; every amplifier has noise and the lower 
bound of possible noises is thermal noise, without which 
the second law of thermodynamics would be overthrown. 
Therefore it appears to be of both theoretical and practi- 
cal interest to calculate the receptivity of boundary-layer 
instability modes to thermal noise and compare its effects 
to e^ predictions. So this paper purports to do. 



II. BASIC FORMALISM OF STOCHASTICALLY 
DRIVEN LINEAR SYSTEMS 

We briefly summarize the theory of stochastically 
driven linear systems so to establish a notation for the 
required formulae. For a forced finite-dimensional linear 
system with state vector f ruled by the evolution equa- 
tion 
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we can define eigenvalues a^''\ and direct u*^*) and adjoint 
v^*) modes (respectively right and left eigenvectors of the 
evolution matrix A). The Laplace transform of eq.(IT]) 
allows us to write its solution as 
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where the expansion of a matrix (with no multiple eigen- 
values) in terms of its eigenvectors (normalized so that 
v(') . u^') = 1) has been used. In the time domain, the 
leading eigenvalue cr'-^-' (the one with the largest real part) 
will dominate for large enough time. Thus, for sufficiently 
large t, 
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If now g(t) is a statistically steady white noise with au- 



tocorrelation matrix (g*(f)g(i')) 
shall have 



G5{t ~ t') [IJ, we 



and the Einstein summation convention is implied on re- 
peated mute indices. Sij are the components of a ran- 
dom white-noise stress tensor whose correlation function 
Sijhk = (sy (i,x)s/,.fe(i',x')) was given in [l| as 

S.jhk = 2iikBT {S.hSjk + 6,kS,h) Sit - t')S{^ - x') (6) 

with kB denoting the Boltzmann constant and T the (un- 
perturbed) temperature. 

As is customary in parallel-flow stability theory, eqs.® 
can be Fourier-transformed with respect to time and 
to the homogeneous directions xi and x^, and thus re- 
duced to ordinary differential equations in the only in- 
dependent variable X2 which we shall henceforth also 
denote as y. We shall follow the convention Ui ^ 
exp (iujt — iaxi —ifixs). Therefore in the following the 
derivative d/dt is understood as a multiplication by iuj, 
the derivative d/dxi as a multiplication by —ia and the 
derivative d/dx^ by — i/3, although where convenient for 
tensorial notation we shall continue to denote them as 
derivatives. 
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Equation ^ is the main result we shall want to extend to 
fluid-dynamic equations in the following. We note that 
in order to determine the mean square value of f at time 
t there is no need to specify the probability distribution 
of g, although in the case of thermal noise this generally 
is gaussian. 



III. FLUID-DYNAMIC EQUATIONS 

We can now apply a similar formalism to the inflnite- 
dimensional system governed by the partial differential 
equations of fluid dynamics. We shall limit our atten- 
tion to incompressible fluids with constant properties in 
a parallel flow. The linearized Navier-Stokes equations 
with an added random stress term can be written as 



duj 
dxi 



= 0, 
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Here ui — U{x2) is a parallel base flow directed along 
the xi axis and depending on the X2 (wall-normal) co- 
ordinate only, Ui denote perturbation velocity compo- 
nents and p the perturbation pressure, p and /i the con- 
stant density and viscosity of the incompressible fluid. 



IV. ADJOINT EQUATIONS 

The homogeneous version of eqs.(l5]), in which they are 
most frequently encountered, is an eigenvalue problem. 
Nontrivial solutions are found when either oj (temporal 
stability problem) or a (spatial stability problem) is al- 
lowed complex values, i.e. more rigorously when either 
the time- or the xi- Fourier transform is replaced by a 
Laplace transform. To study the inhomogeneous prob- 
lem ([SJ, on the other hand, we can obtain a substantial 
simplification by the use of adjoint equations. The ad- 
joint problem is introduced by premultiplying the homo- 
geneous form of ea. l|5ap by a new variable (continu- 
ity adjoint) and the homogeneous form of eg. ((5b)) by a 
new variable (momentum adjoint) ^IS^ and integrating 
over the entire y-domain to obtain the so called Lagrange 
identity: 
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dxi 
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dp 
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dxjdxi 



dxi 



d2/ = 0. (7) 



The Lagrange identity must be satisfied for whatever 
r.h.s. of eqs.®, i.e. for whatever ui and p that obey the 
boundary conditions but not the (homogeneous) direct 
equations. This is imposed by transforming the integral 
through integration by parts until all derivatives of Ui 
and p disappear; the coefficients of these variables in the 
integral are then set to zero, and there results an eigen- 
value problem for the unknowns and qf : 
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Integration by parts at the same time provides the ap- 
propriate boundary conditions for these equations. 

Here we consider the temporal stabiUty problem (sim- 
pler than, but analogous to, the practically more inter- 
esting spatial stabihty problem). If eqs.® are solved for 
a specific (generally complex) eigenvalue u}^^\ whereas in 
eqs.® d/dt is replaced by the Laplace variable a, pre- 

multiplying eqs.® by respectively c^^^^ and q'^^^^ and 
subtracting ([7]) gives 

(a-ic.W)p J ql(%dy = JqlW^dy^ 

^..,dy, (9) 



the last step being once again an integration by parts. 
Since direct and adjoint modes are mutually orthonormal 
with respect to the scalar product / q'j'Ujdy, the l.h.s. of 

eq.® contains the complex amplitude C^^^ of the mode 
u , of frequency in the modal representation of uf. 



explicitly integrated in time. The final result is 



where Rij = —dq^^'^^/dxi. In the time domain 



RijSijdy, 



^ p-i / Kt-r) / R^^s,,Ay dr. (10) 



V. STOCHASTIC DISTURBANCES 

The formalism of the previous sections provides the re- 
sponse of a single mode of our shear flow to a stress forc- 
ing of arbitrary, deterministic or random, nature. Armed 
with these results we can now go back to analyse the ran- 
dom stresses that arise from thermal fluctuations. Since 
Sij is a random white noise, of which only the correlation 
function is known, C'-^-' will be a stochastic variable too. 
The mean square of its value is given by 
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X Rl^{y)Rhk{y')S,jhkdyATdy' dr' ^ ^ . (11) 

On account of eq. ^ two integrals are eliminated by cor- 
responding (5- functions, whereas the exponential can be 
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VI. DISCUSSION AND 
ORDER-OF-MAGNITUDE CONSIDERATIONS 

When the mode with eigenvalue w^^-* is stable 
( Im(ijj(^)) > 0), the r.h.s. of eg. ([72]) converges to a steady 
state for t — > oo, and the fluctuations induced by thermal 
noise have a finite limiting amplitude. When this mode 
is unstable, on the other hand, the perturbation diverges 
in time and no statistically steady state is attained; if 
the instability has a well-defined origin in time (thermal 
noise cannot be switched on and off, of course, but the 
base flow can), ea. p2)) allows us to determine the du- 
ration of the transient until the perturbation attains an 
amplitude comparable to the base flow and transition to 
turbulence occurs. Even more interestingly, an analogous 
process occurs in space: a spatial-stability analysis (the 
details of which will be given in a separate paper) allows 
us to determine the length of the region where laminar 
flow prevails until the perturbation attains an amplitude 
comparable to the base flow and transition to turbulence 
occurs. 

Before proceeding with a numerical calculation, a di- 
mensional analysis may be useful. The dimensions of 
q^^^'' are those of (velocityxlength)^^, i.e. TL~^ if we 
use T and L to denote time and length measuring units. 
The dimensions of Rij are thus TL~'^, and since those 
of {fj,kBT)/p'^ are T^^L"^, the r.h.s. of eq.([12]) is dimen- 
sionless as it should be. If all lengths are normalized 
with a reference length h (a representative breadth of the 
shear flow, such as a boundary-layer thickness or chan- 
nel height) and all velocities with a reference velocity V 
(representative of the base velocity profile, such as its 
maximum value), viscosity becomes encapsulated in the 
Reynolds number Re ~ hV/v (where as usual v — /i/p); 
the coefficient [iiksT) / gives rise to a new dimension- 
less number that we shall conveniently define as the ratio 
of two lengths A//i, where: 



A — ksT/ pv^ 



(13) 



The so defined characteristic length A is a material prop- 
erty; for standard air at 300 K and atmospheric pressure, 
for instance, A = 1.508 x 10"^^ m. Equation (fT^ may 
thus be rewritten in dimensionless form as 
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FIG. 1: Longitudinal-velocity fluctuation iirms vs. local 
Reynolds number, as obtained from a numerical integration 
of the spatial receptivity problem in a Blasius boundary layer 
(solid), compared with the amplification envelope of the stan- 
dard method (dashed). Wrms attains 1% of the freestream 
velocity in the same region where A'^ ~ 9 10. 

where all quantities inside the integral are nondimension- 
alized with the reference length h and velocity V. 

An obvious use of eq.((T4]) is to make a transition pre- 
diction by comparing the value of Crni to the level of 
disturbances that is experimentally observed to occur 
right before the breakdown to turbulence. An order-of- 
magnitude estimate of the expected threshold can be ob- 
tained through the following considerations: the denom- 
inator Im(cL;(^^) is a small number almost comparable 
with Re^^ (the exponent is not quite —1, according to 



the asymptotic theory of the Orr-Sommerfeld equation, 
but close enough for the present estimate) ; the threshold 
r.m.s. velocity perturbation is of the order of 1 2% Q 
while the order of magnitude of Rij is typically larger 
than unity because of the non-normality [lo| of fluid- 
dynamic stability problems. If we assume that all to- 
gether these factors (squared) compensate the remain- 
ing factor of i?e~^, we can roughly equal the exponen- 
tial e~^^™('^*^')*, which would be e^^ according to the 
method, to h/X. For h — 1 mm (corresponding to 
Sqq w 5 mm or L « 1.9 m) in air, we thus obtain 

N = 0.51og(VA) = 9.00. 

Whereas this simple estimate cannot, of course, be 
trusted to three decimal digits, its value is strikingly on 
top of the empirical value of 9 -^ 10 neper used in aerody- 
namics. This is how we got to realize that receptivity to 
microscopic fluctuations could have a significant role in 
causing transition to turbulence. Owing to the presence 
of the logarithm, changing thickness {h) or type of fluid 
(A) just moderately affects this estimate. 

A precise numerical evaluation of the correspondent of 
ea. (|14p for the quasi-parallel case of engineering inter- 
est where amplification occurs in space rather than in 
time is the subject of a separate, longer paper. Prelimi- 
nary numerical results in Figure [1] quantitatively confirm 
that even in the absence of remnant atmospheric sound 
or vorticity (considered until today to be indispensable), 
transition to turbulence on airplane wings would occur 
where it is practically observed to occur. For, microscopic 
fiuctuations of thermal origin set a finite and impassable 
lower threshold to the input noise of this like any other 
amplifier. 
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